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Experimentally the temperature in a Bose-Einstein condensate is always deduced resorting to the 
comparison between the Maxwell-Boltzmann velocity distribution function and the density profile 
in momentum space. Though a successful method it is an approximation, since it also implies the 
use of classical statistical mechanics at temperatures close to the condensation temperature where 
quantal effects play a relevant role. The present work puts forward a new method in which we use 
an ultra-intense light pulse and a nonlinear optical material as detectors for differences in times- 
of-fiight. This experimental value shall be compared against the result here calculated, using the 
Bose-Einstein distribution function, which is a temperature-dependent variable, and in this way 
the temperature of the condensate is obtained. 
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The emergence of the concept of Bose-Einstein con- 
densate (BEG) dates back to the year of 1924 with a pa- 
per by S. N. Bose in which the statistics of the quanta of 
light was analyzed A year later A. Einstein with the 
idea of Bose predicted the occurrence of a phase tran- 
sition in a gas of non-interacting atoms 2]. It can be 
stated that a BEG is a state of bosonic matter, confined 
by an external potential, and cooled to temperatures very 
close to the absolute zero . Under these conditions the 
atoms begin to collapse into the lowest energy wave func- 
tion associated to the corresponding confinement poten- 
tial 0]. 

Since our proposal involves the use of the features of 
nonlinear optics let us explain, very briefly, the main idea 
behind this topic, which emerged in 1961 Q. As a light 
beam propagates through a piece of material the polar- 
ization of its atoms emerges as a response to the excita- 
tion defined by the corresponding electric field. As long 
as the intensity of the electric field remains smaller than a 
certain threshold, the interaction is linear and the trans- 
mitted light has the same frequency as the incident one 
. But if the intensity of the light is beyond this thresh- 
old, then the transmitted light has a spectrum with more 
than one component, among which the double frequency 
is often the most important [6|. 

The experimental achievement of BEG has spurred the 
appearance of many questions in the theoretical realm 
[7| as well as in the experimental level. One of them 
is related to the procedure by which the temperature 
is deduced. Indeed, the temperature is always obtained 
resorting to an absorption imaging method In this 
procedure the corresponding atom cloud is imaged ei- 
ther while it was trapped or following a switch-off of the 
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used trap. The time-of-flight images of the atoms are 
compared against a Maxwell-Boltzmann velocity distri- 
bution function jj], and by this way the temperature is 
deduced. It is noteworthy to mention that in the context 
of this issue the deduction of the temperature is always 
done resorting to a statistics, Maxwell-Boltzmann, which 
is valid in the regime of high temperatures, i.e., far away 
from the condensation temperature [JL It can be argued 
that it is a good approximation (see [i] pp. 27 and [l3| 
pp. 469). Nevertheless, from a conceptual point of view 
this procedure has a severe shortcoming since it resorts 
to a model (Maxwell-Boltzmann statistics) in a region 
in which it loses its validity. This last remark raises the 
following question: could the temperature of a BEG be 
deduced without the use of Maxwell-Boltzmann statis- 
tics? 

The present work provides an affirmative answer to this 
question. Our proposal involves an ultra-intense light 
pulse. This pulse will be divided into two beams, one of 
them will travel through the condensate while the other 
one will not (see Fig. (??)). The difference in time-of- 
flight between these two beams will be detected resorting 
to a nonlinear optical material. Additionally, we deduce 
the expression for the difference in time-of-flight as a 
function of the parameters of the trap, and hence we 
obtain an expression for the temperature. It has to be 
stressed that in this process classical statistics is not used. 
Since the required laser fields generate ultrashort light 
pulses we may detect phenomena involving times scales 
in the range of femtoseconds [6]. This means that our 
procedure offers us the possibility of detecting differences 
in times-of-flight around this order of magnitude. 

Gonsider a laser light beam, with wavelength A, and 
divide it into two parts. One of them will be used as a 
reference beam, and the second one will pass through the 
corresponding condensate. Afterwards, both beams are 
brought together and then the difference in time-of-flight 
between them is calculated. 

Let us now proceed to obtain the time-of-flight as 
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function of the variables of the condensate. The den- 
sity of the gas n(r) can be related to the refraction index 
h{r) and the particle size R [lH as follows (a condition 
to be satisfied is \^Rn{r) < 1) 



n{r) = 1 + \^Rn{r). 



(1) 



The density, which is connected to the eigenfunctions 
of the trapped particles 4'u{'r) and the mean occupation 
number (n^) reads 
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For our case, a Bose-Einstein situation, the occupation 
number is a function of the energy eigenvalues e^, the 
temperature, and the chemical potential ^jl i.e., 
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These last three expressions show that the refraction 
index will depend upon the temperature, and, in conse- 
quence, the difference in time-of-flight will be a function 
of T. This is the point to be exploited in the present work. 
We now consider the case of a Bose-Einstein gas of non- 
interacting and non-relativistic particles immersed in a 
Newtonian homogeneous gravitational field. The inclu- 
sion of gravity has been done because, at least for alkali 
atoms, it has non-negligible effects [H, [l^. Addition- 
ally, an anisotropic harmonic-oscillator confining poten- 
tial is present. At this point we must add that in the 
usual experimental devices the atom clouds are confined 
with the help of laser trapping or magnetic traps [l^l . In 
the case of alkali atoms some of the available confining 
potentials can be approximated by a three-dimensional 
harmonic oscillator (see [l^] pp. 466). In other words, 
our potential represents a very good approximation to 
some experimental cases. 

The frequencies of our harmonic oscillators along the 
coordinate-axes will be denoted by w;, where I — x,y, z. 
In addition, g represents the acceleration of gravity. Un- 
der these conditions the complete potential reads 
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In general the frequencies are not equal (here we as- 
sume uJz l£ ^ ^x), while m denotes the mass of the 
corresponding particles. 

The energy eigenvalues are 
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The density of states is calculated in the usual way 
0i 1 namely, we resort to the continuum approximation 
and consider the number of states proportional to the 



volume of the surface in the first octant bounded by plane 

E=^ ljIx'^ +ujyy'^ +uj'^z'^ +mgz. The derivative, with 

respect to E, provides the density of states (number of 
states per energy unit). 
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The density of the gas (number of particles per volume 
unit) is given by ([2]). At this point we face the following 
question: near the condensation temperature, Tc, how 
many excited state wave functions shall be considered in 
(12)? As a rough approximation we will consider first only 
one excited state. In addition, according to the bosonic 
statistics, below the condensation temperature the av- 
erage number of particles in the ground state reads (for 
the case in which the confining potential is an anisotropic 
harmonic oscillator Q) 
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Additionally, the average number of particles in the 
first excited state is 



Therefore, under these conditions we have that 
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In these last expressions N denotes the number of par- 
ticles comprising the atom cloud and (/)(o) (r) is the ground 
state wave function given as follows 
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The relation between the refraction index h{r) and the 
density n(r) is provided by ([T]) (see [ll| pp. 26) and it 
takes the form 
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The arrival at the nonlinear material of our beams will 
show a difference in time due to the fact that the reference 
beam has always a speed equal to c, whereas the second 
beam will have a lower speed, v, during its pass through 
the condensate. The refraction index is defined as n{r) = 
c/v The velocity inside the condensate is position- 
dependent, i.e.. 
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For the sake of clarity let us assume that the second 
beam during its movement inside the condensate has co- 
ordinates z and y constant. The calculation of the dif- 
ference in optical path requires the knowledge of the dis- 
tance that the light beam travels inside the condensate. 
At this point we introduce a distance parameter stat- 
ing that the size of the condensate along the x-axis is 
fixed by the harmonic oscillator length (see [l^l PP- 467) 
l^. = h/mbJx- 

If Xi and Xf are the coordinates defining this width, 
then 2lx = Xf—Xi. Let us denote by Ai the time required 
by the second beam to move from Xi to Xf. Then 
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Expression is the main result of the present work. 
Notice that it contains a temperature-dependence of the 
time-of~flight inside the condensate. This is the point to 
be exploited. It is noteworthy to comment that Maxwell- 
Boltzmann statistics has not been used at all. 

Our experimental proposal is the following one: A 
ultra-intense light pulse is emitted and divided into two 
parts. The upper beam travels through a rubidium con- 
densate with a certain value of z, say zi, we also keep its 
y-coordinate constant, i.e., y = yi. We choose rubidium 
since its refraction index has already been determined 
by interferometric methods, i.e., the technology that the 
present proposal requires is nowadays available (l5| . Af- 
terwards, both beams are brought together and they im- 
pinge upon a nonlinear optical material. The difference 
in optical lengths between the paths of our two beams 
entails that in general they will not be in phase upon 
their arrival at the nonlinear medium. The measurement 
of temperature will be done using the characteristics be- 
hind the appearance of nonlinear response. Indeed, the 



generation of high harmonics in a nonlinear material de- 
pends upon several conditions, and one of them is related 
to a threshold in the intensity (denoted here by /q) of 
the corresponding light pulse . If the beams are out of 
phase, then the total intensity of the combined beam will 
be smaller than /q and, in consequence, the difference 
in time-of-flight of the two beams implies the absence of 
high harmonics in the response of the nonlinear material. 
If we recover the nonlinear response, then we may assert 
that the beams are in phase. To achieve this we now 
modify the length of the upper arm of the interferometer 
until high harmonics emerge. In this sense we know that 
the difference in optical length is determined by two vari- 
ables: (i) the difference in length between the upper and 
lower arms; (ii) presence of the condensate. With the 
presence of high harmonics we may state that the differ- 
ence in optical length equals a multiple of the wavelength 
of the corresponding light. This allows us to determine, 
experimentally, the difference in optical length induced 
by the condensate. The precision is provided only by 
nonlinear optics, whereas in the usual model it involves 
not only the precision associated to the imaging method 
P but also the error introduced by the use of classical 
statistics. 

We now change the length of the upper arm in such a 
way that at an increase of it equal to A^ the nonlinearity 
appears once again. This last remark implies [16] 
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Quantities A/ and A are detected experimentally, while 
with the help of we determine At. This provides the 
temperature 
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Where I = X — 2lx — AL Let us now estimate the 
feasibility of the idea. Assuming N/L^ ^ 10^^cm~^ 
(see [93 pp. 5), with ^ 10~^cm [ITj], considering a 
neodymium laser A ~ 10~^cto [iB|, and with R equal to 



Bohr's radius, i.e., i? ~ 10 



we obtain as estimation 



A^ ~ 10 •^CTO. Our argument also proves that the exper- 
imental values fulfill the condition behind ([1]), namely, 
X'^Rn{f) ~ 10-2 < 1 iui. 

We now consider a more realistic situation, namely, 
below Tc one finds that q excited states are populated. 
Then ^ becomes 
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In this last expression Hn{x) denotes the Hermite poly- 
nomials [3 and En = huj^n. We now proceed as before. 
Indeed, with (pTjl we obtain the refraction index (resort- 
ing to H])) and calculate the difference in time-of-flight. 
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Let us now evaluate the order of magnitude of the contri- 
butions stemming from those excited states higher than 
the first one. In order to do this consider T = lO^^K 
and also two excited states. Since H{x)s ~ x'^ (when 
X » 1) then772(i + c,/Z^tj2) ^ {l+g/l.^Ljl)'^'' . Therefore 
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Introducing into (j22p the previous experimental val- 
ues we obtain that the contribution of the second excited 
state is several orders of magnitude smaller than contri- 
bution of the first excited state. This argument can easily 
be generalized to g > 3 excited states and it proves that 
for rubidium at T = 1Q~^K only one excited state suf- 
fices for deducing the required time-of flight. 



Without the use of Maxwell-Boltzmann statistics tem- 
perature could also be detected by means of an interfer- 
ometric procedure, but this method has the shortcoming 
of requiring densities which could imply problems to the 
current technology. This last remark is not the proposal 
contained in Q in which temperature was measured with 
interferometry and using Maxwell-Boltzmann statistics. 

Let us mention that our idea could allow us to put for- 
ward a way in which some postulates behind general rel- 
ativity could be tested [l^, i.e., our expressions depend 
upon g. Remember that, up to now, there is no con- 
sistent quantum gravity theory, and that any proposal 
shedding light upon the validity of the postulates behind 
general relativity at quantal realm could be relevant [1^ . 
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